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Abstract 

We compute the critical coupling constant for the dynamical chiral- 
symmetry breaking in a model of quantum chromo dynamics, solving numer- 
ically the quark self-energy using infrared finite gluon propagators found as 
solutions of the Schwinger-Dyson equation for the gluon, and one gluon prop- 
agator determined in numerical lattice simulations. The gluon mass scale 
screens the force responsible for the chiral breaking, and the transition occurs 
only for a larger critical coupling constant than the one obtained with the 
perturbative propagator. The critical coupling shows a great sensibility to 
the gluon mass scale variation, as well as to the functional form of the gluon 
propagator. 
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1 Introduction 



The idea that quarks obtain effective masses as a result of a dynamical breakdown of 
chiral symmetry (DBCS) has received a great deal of attention in the last years |], §]. 
One of the most common methods used to study the quark mass generation is to 
look for solutions of the Schwinger- Dyson equation for the fermionic propagator. It 
is known that above a certain critical coupling (pt c = g^/An) a nontrivial self-energy 
solution bifurcates away from the trivial one. Numerical evaluation of this critical 
coupling in QCD with three and four flavors gives a c ~ 0(1) |3], ||. 

Parallel to the study of DBCS a lot of effort has also been done to obtain the 
nonperturbative behavior of the gluon propagator 0, §^ |7J , and perhaps one of the 
most interesting results is the one where it is argued that the gluon may have a 
dynamically generated mass ||. The study of the infrared behavior of the gluon 
propagator was also performed numerically on the lattice ||], and more recent nu- 
merical simulation give strong evidence for an infrared finite gluon propagator in the 
Landau gauge ||. It is worth mentioning that from the phenomenological point of 
view, the existence of a "massive gluon" may shed light on several reactions where 
long distance QCD effects can interfere, and examples of the possible consequences 
can be found in the literature, see, for instance, Ref. [jlO, |ll|, 12| . 

Much work has yet to be done about the infrared behavior of the gluon propa- 
gator, but it is clear that its implications have to be tested in all possible problems. 
It is possible that the constraint coming from DBCS, and other phenomenologi- 
cal studies jlUL [II], |Tj| will provide a map of the infrared gluon propagator. Since 
the bifurcation point for DBCS was studied up to now with the perturbative 1/k 2 
gluon propagator, it is natural to ask what is going to happen with the infrared finite 
propagators that have been found through solutions of the gluonic Schwinger-Dyson 
equation or using Monte Carlo methods, and, moreover, to look for the consequences 
of different forms of non-perturbative gluon propagators. It is intuitive that the force 



2 



necessary for condensation is going to be screened if the gluon propagator is infrared 
finite, therefore, the actual critical coupling constant should be larger, and this is 
what we will investigate in this work. 

We will present the Schwinger-Dyson equation of our problem, and first we will 
discuss the critical coupling for the linear approximated problem in the case of a 
bare gluon mass. This will teach us on the general behavior of the critical coupling 
constant as a function of the gluon mass. Secondly, we perform a numerical calcu- 
lation of the full nonlinear equation, for two different gluon propagators resulting 
from solutions of the gluon polarization tensor, and one obtained from numerical 
simulation on the lattice. In the conclusions we discuss the differences in the crit- 
ical coupling for each of the "massive" gluon propagators, arguing that its value 
definitively gives information about the infrared gluon propagator. 



2 Quark propagator Schwinger-Dyson equation 

The Schwinger-Dyson equation for the quark propagator in Minkowski space is 

S-\ P ) =^~i\j ^l,S{q)T v {p, q )g 2 D^{p-q), (1) 



where we write the gluon propagator in the form 

9 ^(^^t|W(-^ + f). (2) 

The propagator has been written in the Landau gauge, which will be used through- 
out our work. In the above equations T u (p, q) is the vertex function, and a(—q 2 /A 2 ) 
is the QCD running coupling constant, for which we know only the ultraviolet be- 
havior, and to solve Eq.(l) we make the same ansatz of Ref. ||, || about its behavior 
for the full momentum scale 

2 2 12tt/(33 - 2n f ) 

<*{-q M ) = - 2 ■ (3) 
ln(r + -ft-) 
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Eq.(3) goes continuously to the perturbative result, and has already been used in 
phenomenological applications. 

To proceed further we also need to introduce an ansatz for the quark-gluon vertex 
r^(p, q), which must satisfy a Slavnov- Taylor identity that, when we neglect ghosts, 
reads 

(p-q)^(p 1 q) = S- 1 (p)-S- 1 (q). (4) 

This identity constrains the longitudinal part of the vertex, and if we write S 1-1 ^) 
in terms of scalar functions 

S-\ V )=A{ V ) P :-B{ V ) ) (5) 

we find the solution [[y| 



r"(p )ff ) = A(p 2 )^ + ^-^([A(p 2 )-A(q 2 )]^-[B(p 2 )-B(q 2 )] 

+ transverse part, (6) 

which is a much better approximation than the use of the bare vertex. Assuming 
that the transverse vertex part vanishes in the Landau gauge we obtain 

D^{p-q)T v (q,p) = D^(p - q)A(q 2 ) luj (7) 

and arrive at the approximate Schwinger- Dyson equation 

- 1] i- B(f) = 4 / ~ (<» 

Going to Euclidean space, we will be working with the following nonlinear coupled 
integral equations for the quark wave-function renormalization and self-energy 

16tt r d 4 Q a{{P-Q) 2 /A 2 ) 



[A(P 2 )-1]P 2 = '-fj 



3 J (2vr) 4 $[(P-Q) 2 ] 

P(P-Q)Q.(P-QY 



x [P.Q + 2 



(P-Q) 2 

A 2 (Q 2 ) 



A*(Q2)Q* + B*(Q2) 



(9) 



B^ = ifi* r d4Q ^ p -Q) 2 / A2 ) mq 2 )b(q 2 ) 

{ ) J (2tt)4 $[(P-Q) 2 ] A 2 (Q 2 )Q 2 + B 2 (Q 2 Y l iU J 

where Q 2 = —q 2 and P 2 = — p 2 , and we introduced a function $[(P — Q) 2 ] which, 
in the case of the perturbative propagator, is simply $[(P — Q) 2 ] = (P — Q) 2 , for 
the massive bare gluon it will have the form $[(P — Q) 2 ] = (P — Q) 2 + m 2 , and will 
be a more complex expression in the case of a dynamically generated mass. 



3 The linear problem with a massive bare gluon 

Eq.(|) and Eq.([L0|) possess the trivial solution A(P 2 ) = 1 and B(P 2 ) = for 
small values of the coupling constant. We can also see that B(P 2 ) depends on 
P(P 2 ) at first order, whereas A(P 2 ) has a higher order dependence on B(P 2 ). In 
order to examine the possibility that a nontrivial solution, B(P 2 ), branches away 
from the trivial one at a critical coupling, a c , we examine the so-called bifurcation 
equation |14 |. This involves differentiating Eq.([|) and Eq.(^) functionally with 



respect to B and then setting B = 0. Since the equation for v4(P 2 ) depends at 
least quadratically on B(P 2 ) it will be droped at leading order from the bifurcation 
problem, and we substitute v4(P 2 ) by 1 in the bifurcation equation that will come 
out from Eq.fllQp. We will deviate from the standard bifurcation theory proceeding 
as in Ref. [TjJ, and instead of substituting Q 2 + B 2 (Q 2 ) by Q 2 in the denominator of 



Eq. fllOD , we will replace this term by Q 2 + 5B 2 (0) and define the dynamical fermion 
mass (mf) by the normalization condition 

5B(0) = m f . (11) 

We finally arrive at our bifurcation equation 

5B ^ - ^I^h^°Q^ 



a[(P-Q) 2 /A 2 } 
(P - Q) 2 + m 2 



SB(Q 2 ), (12) 
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where we have already assumed a bare massive gluon. 

Our main intention in this section is to verify the gross behavior of the critical 
coupling constant with the existence of an infrared finite gluon propagator, this 
is the reason for having selected a bare massive gluon in Eq.(^). The details of 
a dynamically generated gluon mass will be left for the next section. Eq.(12) is a 



standard Fredholm equation with a positive kernel, and, requiring 5B(P 2 ) to belong 
to L 2 , the spectrum is discrete with a smallest value a c such that we have the trivial 
solution 5B(P 2 ) = for < a < a c , and the nontrivial one if a > a c . 

We can still make some simplifier approximations before estimating a c , making 
the following substitutions 

a((P - Q) 2 /A 2 ) -> 6{P 2 - Q 2 )a(P 2 /A 2 ) + 6(Q 2 - P 2 )a(Q 2 /A 2 ), (13) 

and 

1 1 g(P 2 -Q 2 )+ 1 0(Q 2 -P 2 ), (14) 



(P-Q) 2 + m 2 g P 2 + m 2 g v ^ ' Q 2 + m 2 g 
which is known as the angle approximation, and introduces an error of about 10% 
in the calculation [|J. Defining the variables x = P 2 /m 2 , y = Q 2 /m 2 , I = A 2 /m 2 , 
k = m 2 /m 2 , and f{P 2 ) = 5B(P 2 )/mf, we obtain 

f(x) = - fdyK(x,y)f(y), (15) 

7T J 

where 

K{x,y) = — ■ — 9(x-y) + — ■ —9{y-x). 16 

x + k y + 1 y + k y + 1 

The kernel K is square integrable 

a 2 (x/£)y 2 



K\Y = / dx dy 



o io (x + K) 2 (y + l) 2 



oo roo 



+ Jo dX L dy ( y + n) 2 (y + ir (17) 
therefore Eq . (p~5|) has a nontrivial L 2 solution for a c on a point set. The smallest 
eigenvalue (a c ) for which Eq.(p!5D has a nontrivial square integrable solution, is the 
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i 


K 


a c 


10 4 


1 


0.6971 


10 4 


10 2 


0.9440 


10 4 


10 3 


1.4853 


10 6 


1 


0.5568 


10 6 


10 2 


0.6607 


10 6 


10 4 


0.9489 


10 10 


10 4 


0.6226 


10 10 


10 6 


0.7822 


10 10 


10 8 


1.2278 



Table 1: Critical coupling constant (a c ) as a function of i = A 2 /m 2 , and k = rrig/mj 
for rif = 4. 

first bifurcation of the nonlinear equation, and satisfy 

- || K ||=1. (18) 

7T 

Table 1 gives the critical value function of £ and K. 

The values of Table 1 were obtained with nf = 4 but they do not change appre- 
ciably as we change nf. As could already be seen in Eq.(|I7]), if we increase the gluon 
masses we can satisfy Eq. fll~8|) only with larger critical coupling constants, and this is 
what we can expect from the numerical solution of the complete nonlinear problem. 
We stress that the values in Table 1, which are approximate solutions of Eq. (|T8"|) , 
can only give a qualitative idea of the problem, because the actual solution, without 
the many simplifications performed until we arrived at Eq. ([17|) , will obviously differ 
from the results of Table 1. 
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4 The critical coupling for infrared finite propa- 



gators 

In this section we solve Eq.flj^) and Eq . (|i~0"|) numerically without further approxima- 



tions. The numerical code we used is the same of Ref. | I6| , and the criterion to 
determine the critical coupling is the one of Ref. Q. With the perturbative gluon 
propagator we obtain (with rif = 4) 



Or 



0.854, 



(19) 



which is compatible with the calculations of Ref. J3], |J. We will solve the gap 
equations with three different propagators which we discuss in the sequence. 

One of the infrared finite propagators found in the literature was determined by 
Cornwall M 



O 2 + 4m 2 (0 2 ) 
<&(Q 2 ) = D^(Q 2 ) = [Q 2 + m 2 (Q 2 )]bg 2 ln[- ' 



A 2 



(20) 



where m 2 (Q ) is the momentum-dependent dynamical gluon mass 

-12/11 



m l(Q 2 ) 



m. 



In 



A 2 



111 A 2 



(21) 



g 2 ~ 1.5 — 2 is the strong coupling constant, and b = (33 — 2n/)/487r 2 is the leading 
order coefficient of the (3 function of the renormalization group equation. This form 
for the propagator was obtained as a fit to the numerical solution of a gauge invariant 
set of diagrams for the gluonic Schwinger-Dyson equation. 

Another infrared finite gluon propagator has been found by Stingl and collabo- 
rators [Q. Its form agrees with that derived by Zwanziger based on considerations 
related to the Gribov horizon |T7[, and is given by 



$(Q 2 ) = D;\Q 2 ) = Q 2 + ^/Q 2 , 



(22) 



S 



where fi is a scale not determined in Ref. ||. It is interesting to note that the 
Bernard et al. || lattice result for the gluon propagator can be fitted by Eq . (|20|) as 
well as Eq. (|22|) . These propagators, apart from some multiplying constant, approach 
the perturbative gluon propagator in the small mass limit. 

Finally, Marenzoni et al. f| also performed a lattice study of the gluon propa- 
gator in the Landau gauge, obtaining for its infrared behavior the following fit 

$(Q 2 ) = D m l {Q 2 ) = m 2 g + ZQ 2 {Q 2 /K 2 )\ (23) 

where m g , Z and rj are constants determined with the numerical simulation. m g is of 
O(A), Z ~ 0.4 and rj ~ 0.5, what is slightly different from the previous propagators. 
The results of Bernard et al. also show the behavior (Q 2 ) v , but with a smaller value 
for 77. 

With the above propagators we computed the dynamical fermion mass as a 
function of the coupling constant. As in Ref. ||] the results were fitted by a function 

h(a) = P(a - a c ) 7 , (24) 

characteristic of a phase transition phenomena. We have not found large differences 
in the values of the critical coupling as we variated the number of fermions, there- 
fore, the fitting will be presented for nf = 4. In Fig.l we plot — 1/ In 5(0) as a 
function of the coupling constant, for the Cornwall propagator (see Eq.( |20])). The 
curves in Fig.l were obtained for m g = 2A and m g = 2.2A, and as expected from 
the example of the previous section if we increase the gluon mass the critical cou- 
pling also increases. These gluon masses are consistent with the values determined 
phenomenologically in the last work of Ref. . The parameters of Eq. ( P4Q and the 
critical coupling are given by 

P = 1.0785, 7 = 0.2535, a c = 0.8692, (m g = 2.0A); (25) 
P = 0.8424, 7 = 0.2682, a c = 1.4211, (m g = 2.2A). (26) 
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As will become clear in the following, not only the value of the gluon mass scale 
is important to characterize the phase transition, but the precise form of the gluon 
propagator will affect considerably the value of the critical coupling. In this case, 
as well as in the next ones, we verified that for small gluon masses we start having 
dynamically generated quark masses for critical couplings quite close to the value 
obtained with the 1/Q 2 propagator (see Eq.(25)). After some value of the gluon 
mass the critical coupling deviates very fast from the value of Eq.(19). An example 
of this behavior is shown in Eq.(26), where the coupling constant is almost twice the 
value of Eq.(25), although we obtained it increasing the previous gluon mass value 
only by 10%! Finally, the chiral symmetry breaking for this propagator was also 
studied in Ref. |E8[| with different results. The main difference lies in the fact that 



in Ref. |L8[ it is assumed a complete cancellation between the coupling constant of 
the vertex function and the coupling in the denominator of Eq.(^). Therefore, the 
mass generation in Ref. [|l^| does not depend at all on the coupling constant, and it 
is far from clear to us that such cancellation should be performed. 

Using the propagator determined by Stingl and collaborators |J, we obtain the 
curves shown in Fig. 2 for /i 2 = 0.25A 2 and /i 2 = 0.30A 2 , and described by Eq. fl24]) 
with 

P = 0.2482, 7 = 0.4784, a c = 2.9038, (fi 2 = 0.25A 2 ); (27) 
(3 = 0.2946, 7 = 0.3362, a c = 6.4720, (fi 2 = 0.30A 2 ). (28) 

Note that the values for the critical coupling constants are quite large. We rely on 
these numbers based on the continuous growth of the coupling constant from a value 
near the one of Eq.(19) for small gluon masses, to the ones of Eq.(27) and (28) as 
the mass is increased. It is known for several other theories with chiral breaking for 
coupling constants of 0(1), that higher order corrections do not modify the critical 
behavior shown by the ladder approximation |TJ|, and we expect the same to hold 
here. Comparing Fig. 2 to Fig.l we see that the dynamically generated mass is much 
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smaller for this propagator, than with the Cornwall one. Performing the calculation 
for fi « (9(3. OA) we do not obtain a significative signal of chiral symmetry breaking, 
i.e. if there is a dynamical mass it is much smaller than A, and do not satisfy our 
numerical criterion to recognize mass generation H. This result is compatible with 



the one of Ref. [p0 |, where it was verified that the quark condensate is consistent 
with zero above a certain critical value of fi for this same gluon propagator. Here 
we foresee a problem for the Stingl et al. propagator, because as shown by Cudell 
and Nguyen JTTJ] we need fi ~ 0(3. OA) to obtain a correct description of diffractive 



scattering with this propagator. We stress that not only the gluon mass scale is 
important to determine the critical coupling constant, but also the functional form 
of the propagator. 

Fig. 3 contains the critical curve for the lattice propagator (Eq. fl23|)) in the case 
of m g = 0.7A, and with 

p = 0.4588, 7 = 0.2870, a c = 3.9712, (m g = 0.7A). (29) 

The Marenzoni et al. propagator gives a value for the critical coupling constant 
which is intermediate between the other two propagators that we discussed up to 
now. If we increase the gluon mass we will also find a point where the symmetry 
is not broken anymore, however, this will occur for larger masses than the one 
predicted in Ref. || (m g ~ A). Comparing all the results it becomes clear that, for 
masses of the same order, the softer is the propagator in the infrared the larger will 
be the critical coupling for chiral symmetry breaking. 

It is known from early studies of DBCS with the perturbative (1/fc 2 ) gluon prop- 
agator [^l|], that the chiral symmetry is broken when the product of the coupling 
constant (a) with a Casimir eigenvalue (CV), which depends on the fermion repre- 
sentation, is larger than a certain critical value. With the introduction of a gluon 
mass scale this is obviously not the case. We have found that the gluon mass scale 
plays an important role in this mechanism, the larger it is the stronger must be the 
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coupling constant to generate dynamical quark masses. In this way, there is some 
chance that the single-gluon-exchange approximation, as presented here, makes no 
sense to obtain consistent solutions of Schwinger-Dyson equations, except for very 
small gluon masses, or for theories with fermionic representations condensating in a 
channel with large eigenvalues Cp. For large gluon masses and fundamental repre- 
sentation quarks the critical coupling becomes quite large, and it is more likely that 
we should consider some effect due to confinement, which could be responsible for 
a correct treatment of DBCS. 



In a related work Papavassiliou and Cornwall |22| also considered the gap equa- 
tion with massive gluons coupled to the vertex equation. They found that non- 
singular solutions of the vertex equation requires large gluon masses, which erase 
the chiral symmetry breaking! In our calculation we introduced a phenomenological 
coupling constant that freezes at low momentum depending on the parameter r of 
Eq. (H) . This parameter is determined when we find the critical coupling a c for each 
gluon mass. Both works can be related if we choose r oc m^/A 2 . With this choice 
when we increase the gluon mass the coupling decreases, and we may never have a 
coupling constant large enough to generate chiral symmetry breaking, arriving at 



the same inconsistency found in Ref. p2[ . We stress that the determination of the 
critical coupling for chiral symmetry breaking, and the verification of the freezing 
of the coupling constant as predicted in Ref. |22| , are crucial tests for the existence 
of a gluon mass scale, as well as can be a good indicator of the true infrared gluon 
propagator. 



5 Conclusions 

We studied the critical coupling constant for the dynamical chiral-symmetry break- 
ing in a model of quantum chromodynamics, using infrared finite gluon propagators 
found as solutions of the Schwinger-Dyson equation for the gluon, as well as one 
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gluon propagator determined in numerical lattice simulations. We first calculated 
the eigenvalue condition for the linear bifurcation equation of the quark self-energy, 
finding that a bare gluon mass scale screens the force responsible for the chiral break- 
ing, and the transition occurs at a larger critical coupling constant if we increase the 
ratio of the gluon to fermion mass. Secondly, we solved numerically the full quark 
self-energy equation for some infrared finite gluon propagators. The result confirm 
our linear approximation, indicating that as we increase the gluon mass the critical 
coupling constant will be larger. We also verified that the functional form of the 
propagator is also important to characterize the chiral transition. 

With the Cornwall propagator (Eq.Q^OD) and gluon masses of the order that are 
expected phenomenologically, we obtain critical coupling constants not far away 
from the one obtained with the 1/k 2 propagator. For this propagator our calcula- 
tion differs from the one of Ref. |Tj| by the reason explained in Section 4. With 
the Stingl et al. propagator (Eq.(^)) the chiral transition will occur only for quite 
large values of the coupling constant. If the gluon mass scale is of 0(A) the critical 
coupling is one order of magnitude larger than the one obtained with the perturba- 
tive propagator. Unfortunately, a phenomenological study of diffractive scattering 
with the Stingl propagator demand gluon masses of 0(3A), for which there is no 
symmetry breaking! This means that this propagator does not represent the actual 



gluon infrared behavior, or the model of diffractive scattering of Ref. [TlJ must be 
modified. The Marenzoni et al. propagator leads to a picture of the chiral transi- 
tion that is consistent phenomenologically, but with a larger value for the critical 
coupling constant. In general, the softer is the propagator in the infrared the larger 
will be the critical coupling. As discussed at the end of last section, for large gluon 
masses, it is possible that confinement effects cannot be discarded, due to the large 
critical coupling constants involved in these cases. The value of the critical coupling 
constant can be used as a tool to study the infrared behavior of the gluon prop- 
agator, and associated with other phenomenological calculations (like the ones of 
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Ref. [pTO] , p] , [T2|]) may provide a map of the gluon propagator for every momenta 
scale. 
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Figure 1: Results of the evaluation of Eq. ([[(]) as a function of the coupling constant 
with the propagator of Eq.(^). We show some of the calculated points, and the 
curve is the result of the fitting by Eq . (|24]) . The calculation was performed for 
rif = 4, the upper curve is for m 2 g = 4. OA 2 and the lower one is for m 2 g = 4.84A 2 . 
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Figure 2: The same as Fig.l for the propagator of Eq.(p2|), with /i 2 = 0.25A 2 (upper 
curve) and /i 2 = 0.30A 2 (lower curve). 
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